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Abstract. In this paper, we study the cosmological dynamics of scalar fields with 
O(A^) symmetry in general potentials. We compare the phase space of the dynamical 
systems of the quintessence and phantom and give the conditions for the existence 
of various attractors as well as their cosmological implications. We also show that 
the existence of tracking attractor in 0{N) phantom models require the potential 
with r < 1/2, which makes the model with exponential potential possess no tracking 
attractor. 



PACS numbers: 0440-b, 9880Cq, 9880Es 
1. Introduction 

Astronomical measurements from Supernovae ^1 12 El and CMB anisotropy jH El 
independently confirm that about two-thirds of the energy density in our Universe is dark 
energy, which has negative pressure and accelerates the expansion of the Universe. The 
most widely studied models for dark energy are cosmological constant and a time varying 
scalar field with positive or negative kinetic energy evolving in a specific potential, 
referred to as "quintessence" with w > —1 [3 IHl UHl UH II2| and "phantom" with 
w < —1 [131 [m HSl Uni HH HHI respectively. Since current observational constraint on 
the equation of state of the dark energy lies in the range —1.38 < w < —0.82 [Tnj, 
it is still too early to rule out any of the above candidates. Complex scalar field as 
quintessence has been proposed in [201 121] and its generalization to the scalar fields 
with 0{N) symmetry as quintessence and phantom in exponential potential have been 
done in I23j . To study the global property of the cosmological system containing 
dark energy, phase space analysis is proved to be a powerful tool. For real scalar field 
quintessence and phantom models in FRW Universe, the phase space analysis have been 
done in [211123 EEl EZI while the analysis has been extended to the 0(A'') scalar fields 
in exponential potentials in Ref . [221 1231 • exponential potential, it has been show that 
the 0{N) phantom does not admit tracking attractor that plays an important role in 
alleviating the fine-tuning problem. In this paper, we study the phase space structure 
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of the scalar fields with 0{N) symmetry in general potentials and discuss the various 
attractors as well as the conditions for their stability. 

2. Phase space of quintessence with 0{N) symmetry in general potentials 

In this section, we study the phase space of quintessence with 0(A'") symmetry in flat 
FRW cosmological background 

ds^ = dt^ - a\t)dx' (1) 

The corresponding equations of motion as well as Einstein equations are just a simple 
generalization of those for exponential potential as shown in |22j . 

h = -y(p^ + p, + R' + ^) (2) 

R + 2HR - — — + V'{R) = 

where R is the radial component of the scalar fields with 0{N) symmetry and all 
"angular" equations have been integrated out and their contributions to the system is 
reflected by the terms associated with S, the integration constant of "angular" equations 
[22] • P7 and are the energy density and pressure of the baryotropic matter and 
p-y = (7 — 1)P7, where 7 is a constant, < 7 < 2. Here, we denote = 8nG. 

Now, introducing the following variables: x = y = ^^-^p^, z = ^ = 

^ ~ ~ KViR) ' ^ ~ ^^v'^R)^'^ ^ ~ ^oga, the equation systemQ becomes the following 
autonomous system: 

^ = ^x[7(l -x'-y'- z') + 2{x' + z')] - (3x - VQz'i - 

% = Ivbil 'x'-y'- z') + 2{x' + z')] - f-Xxy 
^= - 3z + ^^[7(1 -x'-y'- z') + 2{x' + z')] - VGxz^ 

Now, it is straightforward to analyze the critical points as well as their stability. In table 
1, we list the stable critical points, conditions for their existence and the cosmological 
parameters there (Note that we have omitted those unstable critical points). We include 
a detailed analysis for the tracking attractor and the corresponding eigenvalues of the 
linear perturbation matrix in Appendix A and B. 
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Critical points (x,y,z,^) 


Name 


Conditions 








quintessence attractor 


r = 1, A ^ 0, 


< 37 


1 -l + ¥ 


fSle , /37e(2-7e) n n 
V2A'V 2A2 ''^''^ 

0,1,0,0 


tracking attractor 
de Sitter attractor 


r > i,A ^ 0, 7e 

A = 


< 2 


!^ 7e-l 
1 -1 



Table 1. The critical points and their physical properties for 0{N) quintessence 
model. Here, 7e = 2r-i ■ Considering the requirements 7e < 2 and T > i together, we 
have r > ^^"^2^ for the stability of the attractor. 



3. Phase space of phantom with 0(Ar) symmetry in general potentials 

Phantom with 0{N) symmetry in exponential potential has been proposed in Ref . |23j . 
However, in exponential potential, the tracking attractor does not exist and thus make 
the model need careful fine-tuning. In this section, we will consider the model in a 
general potential and give a condition for the existence of tracking attractor. In flat 
FRW spacetime, the equations of motion as well as Einstein equations are 



^ = -Y(P7+P7-^'-^) (4) 

= -3H{p^+p.,) 

V2 

^^ = y[P7-^(^^ + ^) + m] 

By using the same dimensionless parameters as defined in previous section, we can 
reduce the equation system (jl)) to 




% = \yb{l + x'-y' + z')-2{x' + z-)] 

^ = - 3z + ^^[7(1 + x^ -y^ + z^) - 2{x^ + z^)] - ^^Qxzi 
dj _ 

In ~ 
dX _ 

In ~ 




Vqx^x{v - I] 



By comparing the dynamical system of 0{N) phantom Eqs.(jl} with 0{N) quintessence 
Eqs.(j2), one can find that their differences are merely the sign difference of the terms 
associated with kinetic energy. However, these small changes make their dynamics quite 
different. In table 2, we list the critical points, conditions of existence and stability 
as well as cosmological parameters there for 0{N) phantom in general potentials. We 
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Critical points 1 


^x,y,z,0 


Name 


Conditions 






V 2 A ' V 2A2 

0,1,0,0 


,0,0 
^,0,0 


big rip attractor 
tracking attractor 
de Sitter attractor 


r = 1, A 7^ 0, A2 

r < i,A ^ 0, 7e 

A = 


< 37 

< 2 


1 _i_ A! 

A2 /e 
1 -1 



Table 2. The critical points and their physical properties for 0{N) phantom model. 
Here, je = 2r=T 



need to emphasize that in exponential potential, the 0{N) phantom model does not 
admit tracking attractor solution [231 • This is due the fact that the existence of tracking 
attractor requires that the potentials must have F < 1/2 while the exponential potential 
with r = 1 is excluded. Another important point of 0{N) phantom lies in that the 
existence of stable big rip attractor requires A^ < 87, which therefore impose a lower 
bound to the equation of state w > — 1 — 7. The similar property w < — 1 + 7 also hold 
for the 0(A'') quintessence case but the significance is less because this upper bound 
won't be much physical interest. 

From the above qualitative analysis, we need to stress that the critical point 
structure is unchanged from the single-field when the internal 0{N) symmetry is 
introduced since the critical points correspond to vanishing angular components {z = 0). 
However, the angular component will definitely affect the trajectories approaching the 
critical points which we will explore in the ensuing section. 

4. Numerical analysis of the 0{N) quintessence and phantom 

In previous sections, we have studied the phase space of scalar field with 0(A'^) symmetry 
in a flat FRW cosmological background. The critical points indicate that these stable 
attractor phases corresponding to the vanishing of "angular" components resulted from 
the 0{N) symmetry. Thus, what will be the effects of the 0{N) symmetry? In fact, the 
introduction of the "angular" term will change the evolutionary tracks of the dynamical 
systems, which was not well manifested by the above qualitative analysis. In this section, 
we study their dynamical evolution numerically and compare the results with the cases 
of real scalar field models. Here, we choose the potential as the power law of the field 
V{R) = Vq{Rq/ Ry. We choose g = 4 and q = —5/4 for 0{N) quintessence and phantom 
with corresponding F = 5/4 and F = 1/5. The parameters were chosen as Vq = 0.00001, 
Rq = Hq = K = 1 and S = 15. 

From Fig.l to Fig.4, we plot the dynamical evolution of matter, radiation and dark 
energy be it quintessence or phantom with or without 0{N) symmetry. Comparing 
Fig.l and Fig. 2 as well as Fig. 3 and Fig.4, one can find that the internal symmetry will 
affect the dynamical evolution at very early epoch. The early property of quintessence or 
phantom will definitely affect the CMB spectrum and thus provide the possibility that 
how to discriminate them from their counterpart without 0{N) symmetry. Especially, 
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for 0(iV) quintessence, the energy density is not necessarily negligible compared to 
radiation and matter due to the contribution of "angular" term that is proportional 
to a~^, which allows a non-negligible energy density for dark energy at early epoch. 
This is a very interesting feature because if the dark energy is cosmological constant, 
its initial energy density need to be some 100 orders of magnitude or smaller than the 
initial matter-energy density [2H] and thus pose the fine-tuning problem and coincidence 
problem. While in 0{N) quintessence, the "angular" term will contribute a kinetic 
energy that decrease very fast as the universe expands (proportional to a~^). So, the 
initial energy condition could be comparable to the matter which is more natural than 
the real scalar field models. 




Figure 1. Evolution of cosmic parameters for matter, radiation and quintessence 
without O(A^) symmetry. The plot begins with 17m, i = ^r,i = 0.499 and ilfi = 0.002. 
Note also that the label N of horizontal axis has nothing to do with the 0{N) 
symmetry, it is the e-fold number defined as = Ina. 




Figure 2. Evolution of cosmic parameters for matter, radiation and quintessence with 
0{N) symmetry. The tilt of the curve at the early epoch is an important feature of 
0{N) quintessence. The plot begins with flm,i = ^r,i — 0.425 and ilj? = 0.150 
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Figure 3. Evolution of cosmic parameters for matter, radiation and phantom without 
0{N) symmetry. The pfot begins with n^^i = i^r,z = 0.499 and = 0.002 




Figure 4. Evolution of cosmic parameters for matter, radiation and phantom with 
0{N) symmetry. To make the difference between the 0{N) phantom and phantom 
clearer, we alow a negative fin at the initial part, which by definition is positive. In 
fact, in a realistic model, the initial values should be chosen with respect to the energy 
condition p > and the result will be that all curves will go upwards. 



5. Discussion and Conclusion 

In this paper, we study the cosmological dynamics of scalar field with 0{N) symmetry 
in general potentials. We show that although there are merely several sign differences 
between the 0{N) phantom and quintessence system, the dynamical evolution, the 
phase space property and the physical significance are quite different, which lead to 
the difference requirements for the potentials to admit stable tracking attractors. This 
requirement for 0{N) phantom, i.e. F < 1/2 makes the exponential potential with 
r = 1 admit no tracking attractor, which is in agreement with the conclusion obtained 
in Ref.|2Sl- Another important point we need to stress is that in the phantom model, 
the stability of big rip attractor requires that the resulting equation of state w > —1 —7, 
which has very interesting physical implication because this may add a natural priori for 
the lower bound of w and plays a very important role in data analysis of SNelajS]. We 
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also show that for the 0{N) quintessence model, the initial energy density could be larger 
and won't affect the evolution too much because the " angular" kinetic energy decreases 
with a~^. In real scalar field quintessence, tracker mechanism was used to provide a non- 
negligible energy density at early epoch. In the 0(A'') quintessence model, if equipped 
with the same tracker mechanism, it will admit a wider range of initial energy density. 

Appendix A. Analysis for Tracking attractors of 0(Ar) quintessence and 
phantom 

In this appendix, we show the way to analyze the tracking attractors in the dynamical 
systems of 0{N) quintessence and phantom. From the dynamical system listed in 
section 2 and 3, one can not identify tracking attractors. The tracking attractors 
listed in Table 1 and Table 2 are obtained through the following analysis. The tracking 
attractors exist when A is very large and F ~ constant. So, we introduce the following 
dimensionless variables re-scaled by e = j. These variables are X = x/e, Y = y/e, 
Z = z/e and O = ^/e with x, y, z and defined in section 2. Then the dynamical 
system could be rewritten as 



where the plus and minus signs of the last term of the first equation correspond to 
quintessence and phantom respectively. Also, we have dropped the terms containing 
the small quantity e in Eqs. ()A.l|) . Now, it will be straightforward to obtain the critical 
points of the above equation system (jA.lll by setting the righthand sides to zero. The 
results are contained in table 1 and table 2. It is worth noting that the equations for 
0{N) quintessence and for 0{N) phantom are very similar with only a sign difference 
of one term. However, this small change lead to very different expressions for the 
cosmic density parameters at the tracking attractors, i.e. ^ and — ^ for quintessence 
and phantom respectively. These expressions indicate that 7e for 0{N) phantom must 
be negative so that the corresponding makes sense and thus requires that F of 
the potential must be less than 1/2. From the definition of F, it would be clear that 
exponential potential has F = 1 and therefore won't admit tracking attractor in phantom 
dynamical system. 




(A.l) 



dZ 

In 
de_ 

In 



V6X0(F 



1) + - 3Z 
1) 
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Appendix B. The eigenvalues of the linear perturbation matrix for different 
stable critical points 

To determine the stability, we need to linearize the corresponding equation system 
0, © and (jA.lll in a neighborhood of the critical points, and then determine the 
stability of the systems by the eigenvalues of the coefficients matrix. By this way, we 
calculated the various eigenvalues corresponding to each critical points and tabulated 
the results in table Bl and table B2. We need to point out that the eigenvalues 
for tracking attractors of 0(A) quintessence and phantom are the same to each other 
because the linearized equations have only one sign-different term, see Eqs. flA.ljl and this 
difference just does not enter into the corresponding eigenvalues coincidentally. From 
the eigenvalues, we obtain the condition for the stability of the critical points, which is 
contained in Table 1 and Table 2. Note that the eigenvalues for the tracking attractors 
are very complicated, and thus we compute them for the power law potential as we have 
used in section 4. The results are listed in Table B3, from which we know that the 
tracking attractors are stable. 



Critical points (x, y, z, ^) 


eigenvalues 




0, i(A2-6), i(A^-6),A2-37 




of-r -3(2+(2r-3)7e) n , V4+{4-24r)7e + (l+2r)27| ^ 
-^i^^ ^He, 4 [i 1 2+(2r--3)7e J 
3^., ON -3(2+(2r-3)7e)ri V4+(4-24r)7e + (l+2r)272i 
2 We ^J, 4 l-L 2+{2r-3)7e J 


0, 1, 0, 


-3, -3, 0, -37 



Table Bl. The eigenvalues of the critical points for O(A^) quintessence. Note that the 
critical point for the tracking attractor (the second row) as well as its corresponding 
eigenvalues are given in terms of (X, Y, Z, Q) 



Critical points (x, y, z, ^) eigenvalues 



-A 


\/i + f,o,o 


0, 


i(_A2_6), A2-37, |(A2-6) 


V 2 '6' 


^37e(7e-2)^Q^Q 


Q^-p IN., -3(2+(2r-3)7e) n , V4+{4-24r)7e + (l+2r)272i 
■^{^ ^He, 4 [i 1 2+{2r-3)7e J 

ON -3(2+(2r-3)7e)ri V4+(4-24r)7e + (l+2r)272 ^ 
2 we ^j, 4 [-L 2+(2r-3)7e J 




0, 1, 0, 




-3, -3, 0, -37 



Table B2. The eigenvalues of the critical points for 0{N) phantom. Note that, similar 
to the case in Table Bl, the critical point for the tracking attractor (the second row) 
as well as its corresponding eigenvalues are given in terms of (X, Y, Z, 9) 
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Models 


r 


7 


7e 


Re{Xi) 


Re{X2) 


0(iV) quintessence 


5/4 


4/3 


8/9 


-1.1667 


-1.1667 






1 


2/3 


-1.2500 


-1.2500 


0{N) phantom 


1/5 


4/3 


-20/9 


-8.7820 


-2.8847 






1 


-5/3 


-7.2122 


-2.2878 



Table B3. The numerical results of the real part of the eigenvalues for the tracking 
attractors when the background fluid is radiation 7 = 4/3 and matter 7 = 1 in the 
power law potential as we described in section 4. 
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